Equivalence of Classical Skyrmions and Coherent States of Baryons 
I. Constrained Quantization on the SU{2) and 5*0 (3) manifolds 
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In the Skyrme model, the Lagrangian can be quantized in several ways using the collective co- 
ordinate approach. Not all of which produce quantum states that can be interpreted as physical 
particles. For example the SU{2) collective coordinate approach produces both integral and half- 
integral spin and isospin states. Only half of these are the physical baryons. Less well known is 
the fact that it is equally possible to quantize the system using the SO{3) collective coordinates. 
This produces only unphysical integral spin and isospin states. To fulfill the goal of being able to 
express a classical skyrmion as a coherent state of baryons directly in terms of the baryon states, 
surprisingly a combination of both collective coordinate approaches is required. To prepare for the 
subsequent application to skyrmion formation through disoriented chiral condensates in heavy ion 
collisions, the Skyrme model is rigorously quantized using the Dirac prescription for constrained 
systems. This is shown for both the SU{2) as well as the 50(3) collective coordinate approach. 

PACS numbers: 12.39.Dc, ll.30.Rd, 25.75.-q, 03.65.-w 



I. INTRODUCTION 



Four decades ago contrary to the popular trend at the 
time to build theories in terms of fermions, Skyrme suc- 
ceeded in building a model with a Lagrangian density 
based on meson fields. Not that it was particularly dif- 
ficult to construct meson field theories but rather this 
one now known as the Skyrme model has the remarkable 
feature that fermions can come out of boson fields in the 
form of topological solitons |^. Those of the Skyrme 
model arc known as skyrmions. The topological charge 
or winding number B of skyrmions was correctly iden- 
tified by Skyrme as the baryon number and the B — 1 
skyrmion was conjectured to be a fermion |^. It was 
only until much later through the method of Ref. ^ 
that the topological charge was worked out to be equal 
to the baryon number ^. Ref. [|| was also able to 
rewrite the latter in the form of the winding number and 
thus made the connection. As to the conjecture of a 
B = 1 skyrmion is a fermion, this can only be settled 
by considering three or more flavors of quarks when the 
Wess-Zumino effective action [Q was introduced to break 
an unphysical symmetry in the Skyrme model Only 
with this action can one show that under a rotation by 
an angle of 27r, a skyrmion acquires a phase change of 
exp(i7V7r) A B = 1 skyrmion is therefore a fermion 
and a baryon. These distinctive works eliminated any 
doubt on the species of particle that a skyrmion belongs. 
Therefore the Skyrme model in the B = 1 sector became 
a non-trivial model of baryon or nucleon. Many further 
works have been produced using this model and many 
applications were found. For example for studying the 
properties such as masses, magnetic moments, sizes, in- 
teractions, structure, and form factors of nucleons and 
baryons 10|, |ll|, |l^ . Other remarkable points 
about the Skyrme model are that it can explain the fact 



that most spin in a proton is not carried by the valence 
quarks and that there is the existence of the qq sea inside 
the nucleon whereas the more basic non-relativistic quark 
model completely fails in this regard 



16 



There are 



also attempts to go beyond the two-flavor Skyrme model. 
This has proven to be not so simple. There are at least 
two ways to do this. One is to continue with the collective 
coordinate approach initiated in Ref. |^ by embedding 
the flavor SU{2) skyrmion in the higher flavor SU{N) 



theory 17, 1^, the other is to bind a kaon to a fiavor 



SU{2) skyrmion and treating strange quark as distinct 
from the fight quarks |2^, |2) . Unfortunately 

both approaches have their own problems in getting a 
reasonable estimates of the baryon mass spectrum. Some 
reviews of both approaches as well as on other applica- 
tions of the three-flavor skyrmions can be found in Ref. 
p3| , p^ . On two- or more flavor skyrmions these have 
been considered as deuteron or other massive nuclei, one 
can consult for example Ref. [ p5| . 

Our interest in the skyrmion stems from an appar- 
ently rather remote subfleld of physics. While most of 
the works mentioned above were interested in the gen- 
eral properties, static or otherwise, of the baryons, we 
are more concerned with the dynamical production of 
skyrmions in the seemingly rather unlikely context of 
heavy ion collisions. The motivation behind this is the 
possibility of observing skyrmion production. This calls 
for an examination of the skyrmion in a very different 
light. Rather than asking the questions of what are the 
physical properties of the quantum states of the Skyrme 
model, the question we ask ourselves is what exactly is 
the quantum analog of the classical skyrmion solution. 
If one does not know what the classical solution corre- 
sponds to in terms of physical baryon states, there is 
little chance for their observation. 

This possibility of skyrmion production in heavy ion 
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collisions was first raised in Ref. [g6|. Strangely this 
was at least several years before the same possibility was 
considered in jets p7| p8| and reapplied to heavy ion 
collisions (29[ |3^ . Within the context of heavy ion col- 
lisions skyrmions are considered as topological defects 
which might form when the system passes from the chi- 
ral symmetric phase back into the broken phase just as 
strings and monopoles might form during the electroweak 
phase transition in the early universe through the Kibble 
mechanism . However unlike other type of topological 
defects skyrmions are subject to strict boundary condi- 
tions which significantly modify their formation proba- 
bility p2| . This must be kept in mind. Common to both 
heavy ion collisions and the electroweak phase transition 
there is an order parameter or vacuum expectation value 
((/)). Below the phase transition temperature this takes 
on a value and a direction in the normal vacuum. But 
at temperatures above the phase transition, there is suf- 
ficient energy available so that can point in any di- 
rections in the order parameter space. As the system 
cools down different spatial regions can have different 
{4>)'s since they need not be correlated. In heavy ion 
collisions from this point onward there are two scenarios 
which can happen and they are not necessarily mutually 
exclusive. 



The first scenario is the system would try to restore 
the vacuum value of (0), (0) = (cr) for the chiral phase 
transition, by eliminating spatial domains with anoma- 
lous values of {(j)). This is done through pion radiation. 
For example if (f> is pointing in the 7r+ direction in a do- 
main, this is simply eliminated by emitting low energy 
7r+. This has led to the proposal of detecting the chi- 
ral phase transition by detecting the formation of these 
domains. One examines the fluctuation in the yield of 
the three species of pions on a collision event by collision 
event basis. This idea was first raised in Ref. Q 
and gained much momentum through the work of Ref. 
p5| , |3^ . Much of its appeal came from cosmic ray data 
where in a few event there was the absence of 7r°. Such 
isospin breaking occurrences are known as the Centauro 
events. It was thought that one could recreate these Cen- 
tauro events in heavy ion collisions or in high energy 
hadron-hadron collision experiments |3^, 
Soon it was shown that in order to detect these pion radi- 
ations amongst the many copiously produced pions and 
hadrons via other mechanisms, large domains of fixed chi- 
ral orientation pointing away from the vacuum direction 
were necessary |o[]0. Due to the fact that misalign- 
ment of the chiral condensates with the vacuum conden- 
sate is central to this phenomenon, it was coined Dis- 
oriented Chiral Condensates (DCC) soon after the first 
proposal Motivated by the Centauro events, the 

primary observable for DCC is not surprisingly the ratio 
of the yield of charged-to-neutral pions. Much work has 
been done on DCC and especially analysis on how to ex- 
tract information from data. Based on the fact that the 
chiral condensates can point in any directions in the order 



parameter space which is S^, one expect a distribution 



1 dN 



N df 2V7 



(1.1) 



for the fraction of neutral-to-total pion / |39[] . Un- 

fortunately so far all searches for DCC found no trace of 
the existence of these chiral condensates: none in heavy 
ion collisions experiments conducted at the Super Pro- 
ton Synchrotron (SPS) at CERN ||, ||, g and none 
at the Tevatron in Fermilab |Q. Instead of the Eq. (1.1) 



these experiments only found a binomial distribution for 
/ signalling the absence of large chiral domains. 

The second scenario is instead of getting rid of the 
anomalous order parameter by radiating pions away, the 
different chiral orientations of the DCC at the inter- 
faces of the domains may be sufficiently different that 
non-trivial topological defects can develop as the system 
evolves towards the normal vacuum po|| . Such defects 
would be prevented from collapsing altogether and con- 
verting into pion radiations due to the nonlinearity of the 
underlying field theory that dictates the chiral dynam- 
ics. These defects would emerge as anomalous produc- 
tion of baryons and antibaryons via B = 1 skyrmion and 
B = —1 antiskyrmion through DCC. Although baryon 
production through skyrmion is not new, it has only 
been connected to the much discussed DCC and to hy- 
peron data at the SPS recently in Ref. [||, || ||. The 
data at the SPS may be showing an anomalous baryon- 
antibaryon production mechanism which leads to devia- 
tion from the expected rare hyperon yields ^ . 

These works use an approach based much on intuition 
and rough estimates while evidently a more solid theo- 
retical foundation is preferred. For example one does not 
know exactly what the classical skyrmion solutions cor- 
respond to in terms of physical baryons. Therefore one 
has no concrete way to determine what baryons will be 
produced if skyrmions are formed in heavy ion collisions. 
For instance which baryons exactly among the octet or 
decuplet baryons would be produced? If more than one is 
produced, what is the probability of a A being produced 
rather than a E? Why should a proton be produced but 
not a S? Our attempt to concretize the idea of skyrmion 
production and better connecting it to experimental data 
has met with a serious challenge. 

A survey of the existing literature revealed that there 
are various past attempts at identifying the skyrmion or 
chiral soliton with a superposition of baryon states. In 
Ref. H^, ^ m, Q soliton solutions were found from 
the linear sigma model and they were identified with the 
hedgehog baryons. In these papers each hedgehog is a 
three valence quark state surrounded by a coherent pion 
cloud. The known baryon states are then projected out 
of the hedgehog. This method enabled these authors to 
study the static properties of the baryons. It also gave 
the probabilities of the distribution of baryons inside the 
chiral soliton. However this method is based on the linear 
sigma model and cannot be used in the Skyrme model 
because the latter is is essentially the non-linear sigma 
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model with additional terms. Another approach is that 
of Ref. (55[ ^ which identifies the skyrmions as J7(4) 
group theoretic coherent states of the Perelomov type 
|37| , pq] . The group algebra is realized by some auxil- 
iary boson fields and the coherent states are states cre- 
ated by these bosons from the vacuum. To study the 
baryons, one projects out the states with the appropri- 
ate quantum numbers from these bosonic coherent states. 
This method is flexible in that the boson number can be 
equated with the number of colors. Thus one can go 
to the large limit. In spite of some very nice features 
of these approaches, we will adopt another. One that is 
more closely connected to the original Skyrme model and 
the baryon wavefunctions derived from it. The collective 
quantization of Ref. produced a set of baryon wave- 
functions that exist on the compact manifold of . As 
we will explain in this paper, we are of the opinion that 
it is better to identify a skyrmion directly as a coherent 
state of baryons. There will be no auxiliary boson unlike 
Ref. |5^, no explicit coherence in a pion cloud but 
only one directly of baryon states. To solve this problem 
coherent states of baryons on compact manifolds will be 
required. Very fortunately the long duration between the 
skyrmion was first discovered and our arrival at this prob- 
lem has permitted much advances in the field of mathe- 
matical physics. Because of which it is our opinion that 
the time is ripe for tackling this problem using such an 
approach. 

Previously the crucial idea of the classical skyrmion as 
a coherent state of baryons and higher resonances was 
presented in Ref. An outline to the solution was 

given. The main goal is therefore to show how to solve 
this problem in detail. Since this is evidently of inter- 
est in its own right, we will concentrate exclusively on 
obtaining the solution. This will be done in two stages. 
The quantization on the compact manifolds will be in this 
paper. The actual construction of the coherent states of 
baryons will be done in a second paper |Q . Applications 
to heavy ion collisions and other considerations will be 
pursued elsewhere . 

Because of the indices in space and the abstract 
group space in the paper, it is clearest to fix the no- 
tations throughout. We will use /i, v for Lorentz indices, 
6, c,d,e — 0,1,2,3 for indices in the four-dimensional space 
of SU{2), and i,j,k,--- — 1,2,3,... for spatial indices. 
Repeated indices imply the Einstein summation conven- 
tion. 

We will first briefly review the Skyrme Model in Sec. || 
followed by quantization via the collective coordinate ap- 
proach The eigenstates will be obtained which is the flrst 
step in making connection with the physical states. In 
the following sections, quantization using both the SU (2) 
and SO (3) collective coordinates will be shown. While 
the flrst one has been done before, the SO{3) quantiza- 
tion in all its details has not appeared in any literature 
that we know of. Showing them together here permits a 
direct comparison. In both cases, the classical form of the 
theory will be shown flrst, then the Dirac brackets will 



be constructed for a proper quantization of a constrained 
system. The quantum commutators then follow from the 
classical Dirac brackets. Differential representations to 
the operators will be found next and the eigenstates will 
be derived from the Schrodinger equations. Finally in 
the appendices some technical details and wavefunctions 
are given. 



II. THE SKYRME MODEL 

In 1962 Skyrme introduced a Lagrangian density in his 
attempt to unify mesons and baryons. The Lagrangian 
density is 



tT[U''df,U,U''d^Uf (2.1) 



where 



U = exp{iT ■ (t)/U} = {(T + iT- tt)//^ 



(2.2) 



/tt is the pion decay constant and g is now known to 
be the p-n-iT coupling. The flrst term is the usual non- 
linear sigma model and the second term was introduced 
by Skyrme ||l|, It is necessary in order for soliton- 
like classical solutions to the Euler-Lagrange equation to 
exist. To find the solution to the equation of motion, 
Skyrme used the so-called Hedgehog ansatz to impose 
maximal symmetry to the solution. It has the form 



U ~ Us ~ expjiT • r F(r)} 



(2.3) 



where F{r) is a radial function. The energy or mass of 
the static soliton ca n be worked out by substituting this 
ansatz into the Eq. (2.1) 



M = - 



Us) 



+ ^-n^^[2(^)^ 



^]}(2.4) 



In order for the solution to be stable, the energy associ- 
ated with C/g is minimized by a functional variation of F 
with the fixed boundary conditions 



F[r — > oo) 







and 



F{r = 0)=Btt . (2.5) 



This ensures the stability of the solution. The boundary 
conditions are essential for the topological charge to be 
given by the winding number B |Q, ^, [g^, ^ . Because the 
Skyrme model belongs to a class of effective Lagrangians 
that approximate QCD at low energies, the solitons or 
skyrmions of the theory have physical meaning. In the 
introduction B has been identifled as the baryon num- 
ber and a skyrmion is generally considered the "classical 
limit" of a baryon ||l|, g, ||, ^, ^, ||. Unlike a baryon, 
a skyrmion is a classical object without any quantum 
number other than B. Therefore a first step towards 
identifying the classical skyrmion is to quantize it. 
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A. The Hamiltonian 

Quantization of the skyrmion amounts to finding the 
Hamiltonian and the eigenstates. A static skyrmion has 
no time dependence and therefore cannot be quantized. 
Such a dependence must be introduced. One cannot, 
however, have an arbitrary time evolution of U without 
the risk of losing the skyrmion completely or introduc- 
ing unwanted complications during the time evolution 
(though the topological charge is conserved due to the 
fact that the third homotopy group of SU{N) is the 
group of integers ■7T3{SU{N)) — Z). A solution is to 
let the time evolution be restricted to a change from one 
skyrmion solution to another. Given that Us is a solution 
of the Euler-Lagrange equation. 



U's^A{t)UsA^{t) 



(2.6) 



is also a solution for any element A{t) of SU{N) because 
of the form of Eq. {2A). This method of providing a 
time dependence to the skyrmion a nd y et preserving the 
essential form of the solution Eq. (2^) was first shown 
in Ref. [|j and is the well-known collective coordinate 
method of quantization. For simplicity, we will avoid the 
problems of three-flavor skyrmions (l7[ ^ |2^, |2|, 
and deal exclusively with only two flavors. 
To quantize the theory, the Hamiltonian is required 
which can be ob tain ed from the Lagrangian. Inserting 



now U'g into Eq. ( |2.lD and integrating over spatial volume 
gives the Lagrangian 



L = Atr(iit) _ M 



(2.7) 



where A — 

A = 



:5-§VA and 



'df f2sin2^{l + (^)2 + si^} (2.8) 



with f — gji^r. The skyrmion solution Us is now com- 
pletely hidden inside the constants A and the mass M . 
With only the kinetic term depending on the collective 
coordinate A{t)^ the quantization is straightforward. It 
is more convenient to rewrite the element of SU (2) as 



A 



flo 



(2.9) 



Because of the unitary nature of S'C/(2), the components 
af, are subject to the constraint 



ci = + a? = a^ah = 1 



These coordinates are therefore restricted to the surface 
of a four-dimensional unit sphere or S'^ . In terms of a^'s 
Eq. (2.7) becomes 



L = 2Adt,db — Af 



(2.11) 



The conjugate momentum to the coordinate a;, is readily 
found to be TTb — 4Adf,. The Hamiltonian is therefore 



H 



1 

8A 



TTbTTfe -I- M 



(2.12) 



Ordinarily quantization would be straightforward, one 
simply imposes the commutators 



]=0 



[ab,TTc] = iSbc 



(2.13) 



then finds a differential representation for tt^ and solves 
the resulting Schrodinger equation for the eigenstates. 
This would be essentially the approach of which we 
have followed c losel y up to now. The presen ce of the 
constraint Eq. ( 2.1Cl| ) changes everything. Eq. ( 2.13 ) are 
the correct commutation relations only if A specifies co- 
ordinates in Euclidean space R'* which is not the case. 
Because the space is 5^, we will follow the prescription 
of Dirac for quantization in the presence of constraints 

MM- 



B. The Dirac Brackets 

To prepare for the quantization, one first find all the 
secondary constraints from the given number of primary 
ones. Only after all the constraints have been found 
can the Dirac brackets be constructed from the Poisson 
brackets. It is the Dirac brackets and not the Poisson 
brackets that will be carried through into the quantum 
theory as commutators upon quantization Hj. In 
our problem, there is only one primary constraint which 
is Eq. (2.10). To find the secondary constraints, one con- 
structs the total Hamiltonian 



Ht = -^TTbTTft + M + p{abab - 1) 

OA 



(2.14) 



This, in general, is simply the Hamiltonian supplemented 
with the primary constraints via the use of one Lagrange 
multiplier per constraint. To derive the second class con- 
straints, one makes use of the basic Poisson bracket 



{ab, TTc} = Sbc , 

and the recursion formula 

Cj+i = {ci,HT} = , 



(2.15) 



(2.16) 



starting with i = 1, repeatedly until no more new con- 
straints are produced. After dropping non-essential con- 
stant factors, we find 



C2 = hia-bTTb + T^bab) = abiTb = T^bO-b = (2.17) 



(2.10) and 



C3 = -^'^bT^b - ^pabttb = 



(2.18) 



This process terminates at C3 because C2 is reproduced 

at i = 4 



C4 



C2 = 



(2.19) 



According to Ref. ^ Ci and C2 are the second-class 
constraints necessary for constructing the Dirac brack- 
ets while C3 is to provide a solution to the Lagrange 
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multiplier p. To arrive at the Dirac brackets, the non- 
singular anti-symmetric matrix of Poisson brackets be- 
tween second-class constraints must be constructed. This 
is 



m 



{ci,Cl} {C1,C2}\ ^ 2 
{C2,ci} {C2,C2W 1-2 



(2.20) 



Equipped with this matrix, the Dirac bracket between 
any two quantities, say A and B is 

{A, B}d = {A, B} - {A, Ci\mT^^{c^,B} . (2.21) 
m^^ is the inverse matrix of m and is given by 



1/0-1 



2 V 1 



(2.22) 



Applying this result to all combinations of and tTc, one 
gets 

{ab,ac}D = {ab,ac}^0 (2.23a) 

{ab,T^c}D = {Ofc, TTc} - {flfc, C2}m2iHci, TTc} 

= Sbc - abctc (2.23b) 
{7r&,7rc}£> = -{7rfc,Ci}mf2Hc2,7rc} 
-{7rfc,C2}m~^{ci,7rJ 

= —ttbTTc + CLcT^b ■ (2.23c) 

These are the Dirac brackets between the basic positions 
and conjugate momenta on . 



D. Constrained Quantization 

Quantization in the presence of constraints amounts 
to making the change from classical Dirac brackets to 
commutators. That is 



{A,B}i 



-[A,B] 



(2.27) 



The circumflex here over A and B on the right designate 
them to be operators. This reproduces the result for the 
commutators between Sb and tTc in Ref. |Q 

[ab,ac\ = (2.28a) 
[ab,7rc] = i{5bc-abac) (2.28b) 
[7r6,7rc] = i{dcTTb - CLbTTc) ^ iLcb (2.28c) 

Lbc are the operators of the angular momenta. These 
are the basic commutators of the SU{2) collective co- 
ordinate quantization. One can easily find a differential 
representation that satisfies the commutators. In fact the 
representation 



Of, 
TTb 



Qb 



-i{5bc - abac)dc + iaub 



(2.29a) 
(2.29b) 



with any constant a will satisfy them. The choice of a 
and some related issues will be discussed in the Appendix 
0. This is an unfamiliar representation compared to the 
usual derivative form for TTf,. However in spite of this, the 
angular momentum operators are still given by 



Lbc = abTTc - ctcTTb = -i{abdc - acdb) 



(2.30) 



C. Angular Momentum 

On rotation is no different from that in and 
the group is iS'0(4) which has six independent genera- 
tors. These six angular momenta are direct generaliza- 
tion of the three generators of rotation in the familiar 
three-dimensional Euclidean space 



Lbc — ObTTc — acTTfc 

Using C2, one can verify that 



so that the Hamiltonian can also be written as 



^ = -^^bcLbc + M 

Ida 



(2.24) 



(2.25) 



(2.26) 



Based on the geometry of the compact space one can 
expect that the angular momentum is more fundamental 
than the linear ones. This will b e clea rly seen later on. 
This is also the reason why Eq. (2.26) is preferred over 
Eq. (Ol). 



E. Eigenstates 

With the dif feren tial representation of Lf,c, the Hamil- 
tonian in Eq. (^.2(j) becomes 



8A 



- abac)dbdc 
M 



^addd 



M 



(2.31) 



where the can be thought of as the Laplacian on S^. 
Although its form is not of a Laplacian, it does act like 
one on S^. The fact that respects the compact space 
of permits a problem-free (see Appendix ^ differen- 
tiation procedure. As a result, differentiation can now 
be performed in exactly the same way as in Euclidean 
space. The Schrodinger equation can be solved. The 
wavefunctions are polynomials in powers of the combi- 
nations of products of (ofc + iac) Q. For example with 
integral power I 



V'^{ab + iac)' = l{l + '2)iab + iacf 



(2.32) 



where 6 7^ c |^ or with b ^ c ^ d ^ e 

-V^{ab + iac)'-'^{ad + iacf 

= l{l -f 2){ab + iacf-P{ad + iae)^2.33) 
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so the energy only depends on the degree of the poly- 
nomial I and not on the which a^'s that make up the 
polynomial. With the various combinations of {ab + ioc) 
for a given I, there is a high degree of degeneracy in the 
energy eigenstates. To label these states, there must be 
other quantum numbers beside energy. 

While the momen tum operators TTf, acquired unfamiliar 
form in Eq. (2.29b), Ltc remains the same as those in 
Euclidean space. From their combinations spin Ji and 
isospin li operators can be constructed (see Appendix 
^). Using the correspondin g diffe rential operators, it can 
be easily checked with Eq. (2.31) that 







2(J2 



4J2 



(2.34) 



The Hamiltonian then becomes 



H = ^{J^ + P) + M ^ ^J'^ + M ^ ^P + M (2.35) 
4A 2A zA 



and the energy eigenvalues are 



1 



Since 



Ei^—l{l + 2) + M. 



(2.36) 



(2.37) 



the Hamiltonian commutes with spin and isospin opera- 
tors 



(2.38) 



so eigenstates can be labeled by quantum numbers of 
energy the third component of spin m and isospin n 
such as |/, m, n) . For reason that will become clear below, 
it is better to use |Z/2, m,n) instead. 

Acting on a test wavefunction with J3 and , we have 



J'i{ai ± jfls)' 



±|(ai ± ia2) 
±l{ai ± ia2) 



(2.39) 
(2.40) 



J3(ai±m2)' P{aQ±iazY 

= ±^{ai±ia2y-P(aQ±ia^Y (2.41) 

= - p)(ai ± m2)'-^(ao ± la^f . (2.42) 

The degree of the polynomial I not only gives the en- 
ergy quantum number I but also the highest or lowest 
spin and isospin ±1/2 of a (iso)spin multiplet. This can 
also be seen by varying the integer p between and I. 
Also from Eqs. (|3|), ( |2.33| ) and ( |2.34D we see that the 
eigenvalues of and P for the above wavefunctions are 
identical and equal to j{j -I- 1) = l/2{l/2 + 1) or j = 1/2. 
Other examples of the (iso)spin operators actiiig on test 
wavefunctions can be found in the Appendix O Even 
I gives integral and odd I gives half-integral spin and 



isospin states. Both bosons and fermions are therefore 
admitted in this system which was shown to be possi- 
ble long ago in Ref. [^. However not all these states 
are physical. Baryons are fermions and must have half- 
integral (iso)spin. Therefore only odd degree polynomials 
are physical wavefunctions. For example a spin up and 
spin one-half state has I — \ovi^m~n^\/2 and is 
represented by 



1 



(2.43) 



and a state with spin —1/2 and isospin 3/2, I — i ov 
J = n = 3/2 and m = —1/2 has the wavcf miction 

(a|i,i,-i) = [ai + ia2)iao - lasf (2.44) 

TT 

where the Dirac ket notation \a) is used to represent the 
position vector on . Other spin-half wavefunctions can 
be written down. Clearly (ofcitiac)' is the highest or low- 
est (iso)spin state of a multiplet. One can reproduce the 
wavefunctions of the whole multiplet for a given degree 
I by using the (iso)spin raising or lowering operators as 
usual. In Appendix ^ the complete set of normalized spin 
one-half and three-half wavefunctions are given. Quanti- 
zation therefore permits the physical states to be found 
as the eigenstates of the Skyrme Hamiltonian. 



III. SKYRME MODEL WITH SO{3) 
COLLECTIVE COORDINATES 

A. From 517(2) to 30(3) 

The general form of the skyrmion solution introduced 
m Sec. |l| in Eq. (|2^) is invariant under a simultaneous 
SU{2) and spatial rotation. Therefore applying a SU{2) 
rotation alone must generate a rotation in space. This is 
equivalent to mapping fro m 31 1(2) to 50(3). Using the 
trigonometric form of Eq. (2.3) 



Us — 1 cos F(r) + iT ■ f sin F{r) 



(3.1) 



the action of a SU(2) rotation by an element A is 

AUA'' = lcosF(r) +i(AT- rA^) sin F(r) . (3.2) 

1 denotes generally the identity matrix of a group 
throughout this paper. Its dimension depends on the 
context. Here it is the 2x2 identity matrix of SU(2). 
The second term can be further decomposed in terms of 
the Pauli spin matrices as 



At,A^ 



TJRJ^(A) . 



(3.3) 



The coefficients of this decomposition Rij which is a func- 
tion of A form a 3 x 3 matrix which can be considered 
as an element of the S0{3) group. Imagine now that 
the A are the collective coordinates for quantizing the 
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skyrmion. Then the map Eq. (3.3) will permit us to 
rewrite everything in terms of the rotation matrix R. The 
collective coordinates will now have nine components Rij 
(not all independent) instead of the four (also not all 
independent). This method of mapping from SU{2) to 
SO{?)) is known in the literature [pq , p3| . In the previ- 
ous section we have already seen how the SU{2) Skyrme 
theory allows only coherent states formed of both boson 
and fermion states. These states cannot be the equiv- 
alence of a classical skyrmion which must only consist 
of physical baryon states. Let us now construct instead 
the SO{3) equivalence of the classical Skyrme Lagrangian 
and Hamiltonian. Then quantization will give us opera- 
tors and eigenstates in SO {3). 



B. The Classical Theory 



To derive the S0{3) equivalence of Eq. (2.7), we first 
rewrite this Lagrangian using the unitary property of A 
as 



L = -Atr(ylU)2 - M . 
Next we use the identity |Q 

A^A = ir, tr(r,AU) , 
which can be easily verified so that we have 

L = -iA(tr(r,AU)^ 



M . 



(3.4) 



(3.5) 



(3.6) 



To proceed further differentiating both sides of the map 
Eq. (^^) with time, followed by multipl ying A^ from the 
left and A from the right and using Eq. (|3.3|) again to get 



A^An + nA^A = tiRi,{A^)R,,(A) 



(3.7) 



Taking the trace after multiplying on both sides with , 
one arrives at 



Rkj{A^)Rj,{A) = -zefe,ytr(Tj^U) _ 



(3.8) 



The first term in Eq. (3.6) can therefore be replaced to 
give us a Lagrangian in terms of the rotation matrix R 



L = -\\ iv{R-^Rf 



M 



(3.9) 



We have use the fact that Rij{A^) = ^y^(^) here. The 
nine elements of Rij become the (collective) coordinates. 
Similar to the ah these are not independent but are sub- 
ject to the constraints of 5*0(3) 



deti? 

RR-^ 



1 

R-^R = 1 
R^ . 



(3.10a) 
(3.10b) 
(3.10c) 



The superscript T in the last line stands for transpose. 
Because of the last constraint one can get rid of R~^ 



everywhere. More explicitly in component form these 
constraints are 

^ijkRliR2jRzk = ^ijkRilRj2Rk3 — 1 (3.11a) 
^ijkRuRmjRnk ^ijkRilRjrriRkn ^/mn (3.11b) 

RikRjk — RkiRkj = ^ij . (3.11c) 

The first and second constraints are equivalent but the 
second is useful for many algebraic manipulations. 
From the above constraints we have the familiar 

^ik^kj — RkiRkj = —RkiRkj ■ (3.12) 
Therefore the Lagrangian can be rewritten as 

(3.13) 

The conjugate momenta to Rij are 

(3.14) 

The Hamiltonian is readily derived to be 



L — jX RijRij — M . 



— \\Rij . 



H^\X R,jR,j +M^{ HyHy + M . (3.15) 



C. The Dirac Brackets 

As a first step toward quantization, we construct the 
Dirac brackets. To do this, we construct the tota l Hamil- 
tonian in the same way as before from Eq. ( |3.15| ) via the 
Lagrange multiplier method. Let us see how many con- 
straints there are. The first one in Eq. (3.11) appears to 
contain 2x3x3x3 = 54 of them but in fact there is 
only one because of antisymmetry. The first constraint 
is therefore 



ci = eijkRuRijRzk -1 = 0. (3.16) 

The second constraint apparently contains eighteen but 
because of the symmetry and equivalence of the trans- 
posed matrix, there are really only six 



c„ — RikRjk — 5ij — 



(3.17) 



where n = 2, 3, . . . , 7 for {ij) = (11), (12), (22), (13), (23) 
and (33), respectively. That makes seven constraints for 
nine Rij components. But 5*0(3) is a three-parameter 
compact Lie group so there can only be six constraints of 
the type of d for n = 1-7 (constraints that involve only 
Rij and no H^). After dropping cy we have the total 
Hamiltonian 



Ht = \ Hy Hy +M 



(3.18) 



Remember from the SU (2) theory previously that these 
are not the only constraints. The total number of con- 
straints can only be determined by repeated application 
of 



(3.19) 
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for i > 1. This will produce many constraints and we 
are likely to have to calculate the inverse of a very large 
matrix of Poisson brackets in order to construct the Dirac 
brackets. This step is very arduous if not impossible. We 
shall try another strategy. 



1. The Special Case of SO{2): The Constraints 



The simpler C2 and Cg give 

4 = Hn - = 

and 

4 = ni2 + = . 

Continuing with i > 4, we get 



(3.30) 



(3.31) 



To keep the calculation within manageable size, we 
calculate the 5*0(2) Dirac brackets and then extrapolate 
the results to 5*0(3) at the end. The primary constraints 
in this case are 



Cl 



eijRiiR2j — 1 

RikRjk — Sij 



= , 
0, 



(3.20) 
(3.21) 



for n = 2, 3, 4 for (ij) = (11), (12), (22) respectively. Here 
apparently we have four unknowns and four constraint 
equations so that each Rij is fixed and can be calculated. 
However 5*0(2) is a one parameter Lie group. One of the 
constraints must be dependent of the other three and so 
can be dropped. 

There is a choice here of what form of the constraints 
should be used. One could drop one of them, say C4, 
but keeping the first three to ensure that there will not 
be any redundancy, or as is well known there are some 
simple relations between the matrix elements of 50(2) 
which are equivalent to the above constraint equations. 
Let us try the following constraints which are complete 
but simpler 



C2 



i?ll — _R22 
i?12 + i?21 



These evidently imply 



Rii 



p2 

^12 — ^22 



1=0 

= . 



^21=1 



(3.22) 
(3.23) 
(3.24) 

(3.25) 



which shows that they are equivalent to the previous ver- 



The total Hamiltonian can be written as 



Ht = iH„H 



M + Pi{e,jRuR2j ~ 1) 



+P2{Rii - i?22) + P3(i?12 + B21) . (3.26) 

Now using the basic Poisson brackets 

{R,,,Ilki} = S,kSji . (3.27) 

(all other brackets vanish) and after dropping the factor 
Pj/X, repeated application of 



= {c'„Ht} = 



(3.28) 



gives 



2^'y 



(Hiji?2j + R2j^li 
+ Rli^2i + H2,-Rl 



eij{IliiR2j + -RiiH2j) — 



(3.29) 



c'^ = je^jIluU2j-2pi-p2{R22-Rn)+p3{R2i+Ri2) = 

(3.32) 



where Eq. (3.25) has been used. According to Dirac, one 
is permitted to make use of the constraints only after the 
Poisson brackets have been evaluated |Q. Doing this 
now using C2 and C3 reduces Cy to 



-7 — ■J^^ij^li^2j 



Pi =0 



(3.33) 



Next for C5 and Cg , again we impose the constraints after 
the brackets have been evaluated to get 



P2 = 



and 



P3 







(3.34) 



(3.35) 



The constraints C7, Cg and Cg provide solution to the La- 
grange multipliers. There are no first class constraints 
since none of them has vanishing Poisson bracket with 
every other constraints. Therefore the second class con- 
straints are c^, n — 1, . . . , 6. We have a 6 x 6 matrix of 
Poisson brackets of second class constraints {c^, c^}. 



The Special Case of SO{2): The Dirac Brackets 



The Poisson brackets amongst c[, C2 and C3 evidently 
vanish. The same is also true between C5 and Cg. The 
only potentially non- vanishing elements are those Poisson 
brackets involving C4 and the {c2,c'^) and (c'^^c'q) pairs. 
Working them out explicitly, the matrix of Poisson brack- 



ets m,; 



{4: 



c;-} is 




(3.36) 



This can be easily inverted to give 



0-1 On 

0000-10 

00000-1 

100000 

010000 

1 0' 



(3.37) 



Tr ying out a few Dirac brackets using the definition 
Eq. ( 2.21 ) produces vanishing brackets between the com- 
ponents of the coordinates 



{Rij,Rki}D — 



(3.38) 
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Between coordinate and momentum pairs, we can have 
for example 

{i?ii,nii}n = 1 - {_Rii, C4}to4i^{ci, Hii} 
-{i?ii,4}m52i{4,nn} 

— \(\ — R22R22) = 5(1 ~ -Rii^ii) , 

(3.39) 

{i?ll,ni2}D = -{i?ll,C4}m4i^{c'i,ni2} 

= \R22R21 = — ^^ii-Ri2 , (3.40) 

{i?12,ni2}_D = I - {Rl2-,c'^}mll{c[,Iil2} 
-{i?12,cam6-3Hc3,ni2} 

— 5(1 — R21R21) = 5(1 — R12R12) , 

(3.41) 



and also 

{Ri2,nn}D = -{Ri2,c'^}m^^^{c[,Un} 

= \R21R22 = —\R12R11 ■ (3.42) 

The same must hold if the indices have been interchanged 
(1.2) — > (2,1) everywhere. One could now attempt to 
deduce the general expression for the Dirac brackets be- 
tween Rij and 11^;; 

{Rij,Tlki}D = SikSji - ^{6,kSji + R^JRkl) + ■■■ ■ (3.43) 

The first term is from the first P oisson bracket in the def- 
inition of the Dirac bracket Eq. ( ^.2l[ ) and the third term 
seems to hold everywhere. The second term is there so 
that this also agrees with Eqs. ( 3.39| ) and ( ^.41 ). Lastly 
there might be other terms hence the ellipses. 

Other Dirac brackets between i?,, and Hi./ can be de- 



duce using the constraints. From Eq. ( 3.3E ) 
{-Rii, n22}D = {-R22, nn}i3 - i(l - R11R22) , (3.44) 



the first term of which deviates from Eq. dHI). Then 
from Eq. ^Idj ) 

{i?22,ni2}_D = — {^11, n2i}£) = — {-R22, n2l}D 

= —\R22R12 — \R11R21 — \R22R21 ■ 

(3.45) 



These are all c onsistent with Eq. (3.43). Next there is 
also Eq. ( 3.41 ) which can be combined with the con- 
straints to give 

{i?12,n2l}D = {i?21,ni2}D - -i(l +i?12-R2l) (3.46) 

and 

{i?2i, n2i}D = i(l - R21R21) ■ (3.47) 

Of these equations the second is consistent with our guess 
in Eq. (3.43) but not the first. This seems to hint at a 



term of —SuSjk/'^ whe n i ^ k or j I otherwise this 
would contradict Eq. ( |3.39| ). Then Eq. ( ^.44| ) suggests 
another term of the form SijSki/2 under the same condi- 
tion of z ^ fc or j ^ Z. The sum of these two terms would 
ensu re all conditions are met. This requires extending 
Eq. to 

{Rij,Iikt}D = 3ikSji — ^{SikSji — Sij6ki + SiiSjk + RijRki) ■ 

(3.48) 

This is consistent with all the above Dirac brackets. Sim- 
plifying as much as possible finally gives 

{R^j,Tlkl}D = ^{SikSji +eikeji - RijRki) ■ (3.49) 

We have made use of the identity 

fifcej; = S.jSki - SaSjk (3.50) 

which is only a special case of the three dimensional ex- 
pression for eij3emn3- 

There still remain the Dirac brackets between mo- 
menta to be worked out 

{n.y , = -{Ilij,c[}m~l{c'4,Ilki} 

-{U,j,c'^}m^^^{c[,Uki} . (3.51) 

For example between equal-index pairs, all Dirac brackets 
vanish 

{nii,nii}i3 = {n22,n22}D - {nii,n22}D = o . 

(3.52) 

Between equal pairs, they also vanish 



{ni2,ni2}i5 = -{ni2,n2i}i3 = o . 

Next between unequal-index pairs 



(3.53) 



{nii,ni2}z3 — 5(^221121 — -^211122) 

= -i(i?iini2-i?i2nn) . (3.54) 

To arrive at the last expression, one uses c'2 and Cg. This 
last line seems to be more closely related to the LHS. 
From here the rest of the brackets can be deduced 

{nii,ni2}D = {n22,ni2}i5 = -{nn,n2i}z5 

= -{Tl22,n2i}D^0. (3.55) 

Aiming now for a general expression, from these exam- 
ples we can infer that the Dirac brackets between conju- 
gate momenta are 



{U,,,nki}D = -^R^jTlki - RkiTl, 



(3.56) 



This is satisfied by the Dir ac bra ckets above. This can be 
further verified using Eq. ( 3.49 ) and the Jacobi identity. 
The latter is one of the basic requirements for the Dirac 
brackets [Q. The Dirac bracket between conjugate mo- 
menta are in general non- vanishing. This is different from 
the familiar examples when the different components of 
the momenta are independent. 
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3. Extrapolating The Dirac Brackets from S0{2) to 5*0(3) 

Returning now to iS'0(3), we still have the original pri- 
mary constraints of c^, i = 1-6 which ultimately dictate 
what the Dirac brackets will be. We start by considering 
the simplest ones — those involving only components of 
R. The vanishing of the Dirac brackets between them 



{Rij,Rki}D = 



(3.57) 



is expected to continue to hold for general SO{N). This 
can be understood as follows. Since the only non- 
vanishing Poisson brackets between i?y and the con- 
straints are those c„ which involve Iljj, these must all 
be secondary constraints derived from the c,, i < 6 via 
{ci,HT}- In this way, each primary constraint will lead 
to a secondary constraint c„, 6 < n < 13 such as 

C7 = eyfe(niii?2ji?3fe + Rli^2jR3k + -^11-^2^113^) = 

(3.58) 

and 



UikRjk + RikTljk = 



(3.59) 



for 8 < n < 12 with = (11), (12), (13), (22), (23). 
The next level of constraints c„ with n > 12 generated 
from these will therefore be equations involving the La- 
grange multiplier pi, i < 6 since the presence of IT will 
bring the multipliers into the equations. We expect no 
more new constraints after that. The matrix of Poisson 
brackets between second class constraints = {ci,Cj} 
will be a 12 X 12 square matrix. It must have a 6 x 6 
block of zeros in the upper left-hand corner between row 
one to six and column one to six because primary con- 
straints have mutually vanishing Poisson brackets. After 
taking the inverse, this block will reappear on the lower 
right-hand corner of m^^. All cofactors in this 6x6 cor- 
ner vanishes because of the block of zeros in m. But the 
only way to have non-vanishing Dirac brackets between 
components of R is if the entries in this block are not all 
zeros because Cp} only if 7 > p > 12. Thus the 
R components must commute with e ach other. 

Let us now extrapolate Eq. ( |3.49 ). The only part of 
this equation that cannot readily be extended to iS'0(3) 
is the second term with the antisymmetric tensors. In 
5*0(3) the structure constants are the e^fc instead of . 
The latter is implicitly equal to eij3 when i, j a re re- 



stricted only to 1 and 2. The eik^ji in Eq. ( 3.49 ) must 
be replaced by eika^jis- But this would leave us with the 
two 3 indices when there is no reason for 3 to be any 
more special than either 1 and 2. If one remembers that 
5*0(2) is a subgroup of 50(3) which is equal to the for- 
mer when the constraint i?33 = 1 is imposed. With this 
piece of information, one can now remove the 3 indices 
by first identifying Ciksejis with Siks^jiaRss- Then one 
can replace this with eikmf-jinRmn when the first pair of 
indices are no longer restricted to 1 and 2. The Dirac 
brackets become 

{Rij^^kl\D — hi^ikSjl+^ikm^jlnRmn — RijRkl) ■ (3.60) 



This must hold because one could imagine using the in- 
dices, say 2 and 3, for 50(2) instead of 1 and 2 then this 
would reduce to the correct form even for this artificial 
case when Rn = 1 and Ru — Rn — for i — 2,3. 

Lastly the brackets Eq. ( |3.56 ) seem to be ready for 
generalization to 50(3) by letting the indices running 
from 1 to 3. Are these the correct Dirac brackets for 



50(3)? In fact for 50(2), the form of Eq. (|3.56| ) is not 
unique. In this special case, one coul d equ ally rewrite it 
after some experimentation with Eq. (3.54) using c'2 and 
C3 as 

{Ilij,Ilkl}D = jSik{Rml^mj — Rmj^ml) 

+ j5ji{Rkm^i7n — RiTTi^ki) .(3.61) 

Which of these are the correct Dirac brackets? It can be 



easily checked by using the Jacobi identity that Eq. ( |3.6l| ) 
are the correct ones. 



D. Spin and Isospin Operators 



The Lagrangian Eq. (3.9) manifests symmetry under 
the left and right 50(3) rotation. An infinitesimal trans- 
formation with parameter rjk, where k index labels the 
direction of rotation, gives 

Rij — + Rij + VkSk Rij — Rij + Vkf^kiiRij (3.62) 
Rij — > Rij +rikSk Rij = Rij + VkRii^kij ■ (3.63) 

Because R has only time dependence, these generate the 
conserved left and right Noether charges but no three- 
current 

dL 

Qi = TTa — — UjkeijlRlk — SijkRkl^jl 



dRjk 

dL 

dRjk 



(3.64a) 



Rjk — ^jkRjl^ilk — Sijk^lkRlj 



(3.64b) 



Time evolution of the primary constraints in Eq. (3.11) 
leads to a set of secondary constr aints [|64| |6^ . For ex- 
ample the second equation in Eq. (3.11) gives 



Rik^jk + ^ikRjk — Rki^kj + ^kiRkj — . (3.65) 

This together with the primary constraints allow us to 
derive 



QfQf = QfQf = 2n„n, 



. (3.66) 

The Hamiltonian acquires yet other equivalent forms 



2A 



M 



2A 



M 



(3.67) 



The significants of the conserved charges will become 
clear in the next section. 



11 



E. The Connection between the SU{2) and SO(3) 
collective variables 



To get an idea of the physics of the conserved charges, 
we now derive the connections between the SU(2) con- 
jugate pairs (abjTTf,) and those pairs , Ilij) of 50(3). 
From the map Eq. (3^), one can deduce that 



(3.68) 



To make the connection more explicit, one can use the 
component form of A in Eq. (2.9) to get the expression 



2aiaj + Sij{2aQ — 1) — 2eijkaoak 



(3.69) 



This expression together with Eq. ( |3.14 ) give us the con- 
jugate momentum Ily in terms of at, and Trt, 



n 



-<^i]k{aaTrk + T^oak)} ■ (3.70) 



Substituting these functions into Eq. (|3.64 ) and after 
some algebra, we have 



Qf 



^{aiTTo - aoTTi - eijkajTTk) (3.71) 
^{aiTTQ — aoTTi + eijkajTTk) ■ (3.72) 



Com pari ng the expressions for Ji and li in Eqs. (C4) 
and ( |C5[ ) (temporarily replacing tti, with —idb to see this) 
reveals that 



Qi — Ji : Qi — 



(3.73) 



The conserved charges are actually spin and isospin. 
Note that all variables here are classical and commut- 
ing so that there is no ambiguity as to the ordering. For 
this reason, we will refrain from carrying the connections 
here to the quantum level when one does not know which 
operator ordering is the correct one. However one still 
identifies the charge operators with spin and isospin op- 
erators as shown here. Therefore the Hamiltonian can be 
written in terms of spin and isospin as in the SU (2) case 



H 



1 

2A 



.P + M 



1 



M . 



(3.74) 



F. Constrained Quantization 

Quantization is accomplished by simply replacing the 
Dirac brackets by —i times the commutators. There- 
fore 

[R^J,Rkl] = (3.75a) 

[Rij.tikl] — \i{5ikSjl + eikm^jhiRmn — RijRkl) 

(3.75b) 

[nyjllfci] — ji^ikiRmli^mj — Rrnjiirnl) 

+ jiSjl{Rk,niii7n ~ Rimtlki) ■ (3.75c) 



A simple multiplication and differential representation 
of the operators are 

Rij — Rij (3.76) 

liij = —^i{5ik5jl + ^ikm^jlnRmn — Rij Rkl)dkl + iaRij 

(3.77) 

where the shorthand dij = d/dRij is used. Promoting 



Eq. ( |3.64 ) first to operators, then after some algebra us- 
ing the identities 

^mja^nlbRklRab Rjnn^jk Rjn^mk 

(3.78a) 

^mla^nkbRljRab = RmnSjk — RmkSjn (3.78b) 

and the above representation, it can be shown that 

Ji — ■^eijk{Rjitlki — Rkiiiji) 

= -\i^i3k{R]idki ~ Rkidji) (3.79) 

li — \Ujk{Rlj^lk ^ Rlki^lj) 

— ~\'''^i]k{Ri]dik — Rikdij) ■ (3.80) 

One can easily check that they satisfy all the commuta- 
tion relations expected of the spin and isospin operators 
(these are of course the same as those for SU (2) see Ap- 
pendix 

G. Eigenstates 



From Eq. ( p. 74 ) and the J and / operators, the differ- 
ential form of the Hamiltonian is 



2A 



- 2Ri. 
M . 



RimRkjdkm)dij + M 



(3.81) 



To find the wavefunctions, it is more convenient to use 
specific combinations of the components. They are 

(3.82a) 
(3.82b) 
(3.82c) 
(3.82d) 
(3.82e) 
(3.82f) 
(3.82g) 



R++ 


= Rii 


— R22 + 


l{Rl2 


+ R21) 


R+- 


= Rii 


+ R22 — 


l{Rl2 


-R21) 


R-+ 


= Rii 


+ R22 + 


l{Rl2 


-R21) 


R— 


~ Rii 


- R22 — 


i(i?12 


+ R21) 




— Rl3 


±iR23 








— R31 


±iR32 






R°° 


= R33 









Using these, it can be easily verified that 



(3.83) 



where I e Z with the exception of p = (7 = when I > 1, 
or 

- v2(i?±9)'-'-(i?±0)'' = l{l + l)(i?±9)'-'-(i?±0)'- (3.84) 
with r < I. The energy eigenvalues are then 

Ei = :^l{l + l) + M . (3.85) 
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It can thus be established that the wavefunctions are 
polynomials in the products of the components of R. 
They are also eigenstates of the third component of spin 
and isospin 



APPENDIX A: 



MR"")' 

IsiR^'^y-^'iR^'^Y 



pliRP'iy (3.86) 
qliRPiy (3.87) 
±Z(i?±«)'-''(i?±°)'- (3.88) 
g(/-r)(i?±»)'-'-(ii'±°)'' . 

(3.89) 



The polynomial degree I dictates the energy, spin and 
isospin eigenvalues just as in the SU{2) theory. This 
time, however, there are only integer values of I (functions 
with half-integer I are not defined). The Hilbert space 
of this system consists purely of bosonic states |j, to, n) 
where j, m,n G Z. These are not physical states expected 
from the Skyrme model. 



IV. SUMMARY AND OUTLOOK 

In this paper we have achieved the quantization of the 
Skyrme model on the as well as on the SO{3) compact 
manifolds. These have been done in a rigorous manner 
using the Dirac prescription for a constrained system. 
While the results from the SU{2) collective coordinate 
quantization are mostly known, this is not true for the 
same on the 5*0(3) manifold. For example, the commu- 
tators Eq. (3.75) in this case have not appeared in any 
literature that we know of. This may be related to the 
fact the eigenstates are not physical and therefore there 
is little incentive to pursuit this second method in any 
detail. This will turn out not to be completely true. In 
the next paper to come and indeed was mentioned in 
p9| , we will need the SO{3) quantized operators, if not 
the unphysical eigenstates, to construct coherent states 
of baryons. The method of coherent states on a compact 
manifold is a mixture of methods used in the context of 
functional analyses and in the studies of the classical 
limit of quantum gravity |6|, ^ |7|, p, |74[ All 
these will be shown in the work to come 1 16 Oil . 
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CHOICE OF THE SECONDARY 
CONSTRAINT 



Quantizing the theory produces certain ambiguities as- 
sociated with the constraints. These ambiguities have to 
do with the ordering of the classical c-numbers turned 
quantum q- numbers. For example while the primary con- 
straint 



ci = flftflb — 1 = 



(Al) 



does not pose a problem, the secondary constraint 
Eq. (2.17) can be expressed as any of the following: 



C2 = abitb = 



C2 = TTfcafc = 



or 



C2 



TTbab) 



. 



(A2) 



(A3) 



(A4) 



They are not equivalent because of Eq. ( 2.28b ). Let us 
not insist on any one of them. From Eq. (2.28b) they 
imply 



-z(rf-l) , 
i{d~l) , 
iz(d-l). 



(A5a) 

or (A5b) 

T^'bO-b ~ —hi{d — 1) , (A5c) 



respectively. We have decided to introduce the dimension 
d during the intermediate steps of the calculation which 
can be set to d = 4 at any time. We have seen in the 



differential representation of TTf, in Eq. (2.29b) that there 
is a free parameter a. Using this representation and the 
above equations, one can see that a is in fact not really 
a free parameter. Its value is actually restricted by how 
the secondary constraint is implemented. In the three 
cases, we have 



a = , 
a ~ d — 1 , 
a = i(rf — 1) , 



or equivalently 



(A6a) 
(A6b) 
(A6c) 



(A7a) 
(A7b) 
(A7c) 



respectively. Among the three choices, the third one 
seems more natural in the sense that Eq. (A4) is a direct 
consequence quantum mechanically of the time-evolved 
form of Eq. (|A1|). 

In the main text, the readers might have noticed that 
we always expressed the classical Hamiltonian in terms of 
first and then upon quantization substituting every- 
thing with their corresponding operators. This is slightly 



TTfc = 


-i[Sbc 


— abac)dc , 




TTb = 




- abac)dc H 


- i{d - l)ab , 


TTb = 




- abac)dc H 


- ^i{d — l)db 
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different from quantizing with the classical Hamiltonian 
expressed in terms of the Tr^Tr;,. We have been careful to 
avoid the latter because of the non-uniqueness of the dif- 
ferential representation of TTf,. On the other hand the dif- 
ferential representation of Ji is independent on the choice 
of a. We feel that the Hamiltonian should be free from 
such ambiguity and in any case Ji are more fundamental 
on than TTf,. In favoring Ji over tt;,, the Hamiltonian 
operator is 



d+1) 



(A8) 



The last equality is more cumbersome looking. The extra 
terms are needed so that different values of a will give 
the same Hamiltonian because the Ttb will change with a. 



APPENDIX B: QUANTIZATION ON 



Quantization of the Hamiltonian Eq. (2.12) is prob- 
lematic because the coordinates A lives on S'^ which is a 
curved compact space. It is tempting to follow the usual 
procedure of c anon ical quantization and use the commu- 
tators in Eq. ( |2.13| ) 







i5i. 



(Bl) 



This would give the familiar differential representation 
for T^b 



-Kb = -iT^ = ~lOb 

dab 



(B2) 



The problem of following this standard procedure is that 
in the subsequent use of TTf,, it is always necessary to keep 
in mind that the derivative is meant to be on S^. Thus 
the straightforward rules of differentiation does not nec- 
essarily apply and indeed may be modified. For example 
while dbUc = 5bc in I^", this is not true on which can 
be very confusing. In fact one must use the strange result 
dbttc = Sbc — abttc to get the ri ght a nswer. The problem 
is that the commutators Eq. (2.13) do not respect the 
constraint of this system which is abO-b = 1- A correct 
quantization procedure must yield 

[abdb - i, O] = [ubCLb, O] = , (B3) 

where O is any operator, whereas Eq. ( |2.13| ) gives 

[dfeflb, aj = , [abttb, TTc] = 2iac ■ (B4) 

The first one of these is fine but not the second. In the 
presence of constraints, one quantizes the Dirac brack- 
ets instead of the Poisson brackets. In short the Dirac 
bracket is the Poisson bracket but with some subtractions 
so that every constraint of the system is respected by the 
Dirac bracket. Then upon quantization they are auto- 
matically respected by the commutators. In essence this 
is not dissimilar to the way of constructing a vector from 
a given vector so that the result is orthogonal to a third. 
In the main text we showed how to do this properly. 



APPENDIX C: ANGULAR MOMENTUM, SPIN, 
ISOSPIN OPERATORS AND THEIR 
EIGENFUNCTIONS 

In the text we gave the six generators of rotation in 
or those of the group 50(4) 

Lbc = abTTc - flcTTf, ~i[abdc ~ acdb) ■ (CI) 
Rewriting them more conveniently in three-vector form 

1 



Ki — ioi 

they obey the following commutators 

[Li,Lj] ~ iSijkLk 
[Li,Kj] = lEijkKk 
[Ki, Kj] = iSijkLk ■ 



(C2a) 
(C2b) 



(C3a) 
(C3b) 
(C3c) 



From them one can construct generators of two separate 
copies of SU{2) 

Ji = i(Li + k^) = ^i(aidQ - a^d., - Ei-jkajdj^ (C4) 



\{U-ki 



aido - Sijkajdk (.C5) 



They are the spin and isospin operators with the familiar 
commutators 



[J^,^/^] — i^ijk'-^k 
[I-iilj] = iSijklk 
[J,,/,] = 0. 



(C6a) 
(C6b) 
(C6c) 



One can form raising and lowering operators for spin and 
isospin as usual 



J± = Ji± iJ2 , I± — Ii± ih ■ 
Their commutation relations are 

[J+J-] = 2i3 , [i3,i±] = ±i± : 
[/+,/_] = 2/3 , [hJ±]^±i± 



(C7) 



(C8a) 
(C8b) 



which follows naturally from Eq. (C6). As is well known 



the squared operators can be expressed in terms of them 



J' = J-J+ + J3(J3 + 1) = J+J- + J3(J3-1) 



(C9a) 



= /-/++/3(/3 + l) = /+/-+/3(/3-l) 



(C9b) 
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APPENDIX D: WAVEFUNCTIONS ON 

The operators J3 and I3 acting on some general rep- 
resentative wavefunctions is enough to show the range 
and relation of the spin and isospin quantum numbers 
to the polynomial degree I. Acting on wavefunctions of 
(flb + ittc) raised to some power give 



/3 

(a|l, 0,+l) = -i — (ai +m2)(ao -iag) (D17) 
(a|l, 0, 0) = ^{ial + al)-ial+al)}iB18) 

TX\J 1 

(a|l, 0,-1) =-i (ai-m2)(ao + ia3) (D19) 

TT 



J3(ai±m2)' = ±5(ai±ia2) 



Jz{aQ±ia-i) = ±2(00 iiaa) 
I 



/3(ao±M3)' = =F5(ao±ia3)' 
and on a product of two different (a;, + iac)'s give 



(Dl) 
(D2) 
(D3) 
(D4) 



= ±i(ai±za2)'"^(ao±M3r (D5) 
= ±(|-p)(ai±ia2)'-^(aoT*a3f (D6) 

= ±(|-p)(ai±za2)'"^(ao±M3)P (D7) 

/3(ai ±m2)'"P(ao Tias)^ 

= ±|(ai±za2)'"^(aoT*a3)^ (D8) 

where I > p and p > 0. Varying p between and I shows 
that Z/2 is the largest modulus of the third component 
of spin and isospin of the multiplet. 

In flavor SU(2) Adkins et al ^ showed that one could 
obtain the normalized wavefunctions on . We give the 
complete set for j = 0, 1/2, 1 and 3/2 wavefunctions. 
The relative phases of these wavefunctions in a multiplet 
are important. For j ~ this is simply a constant 



(a|0,0,0) = 


1 

7rV2 ■ 




(D9) 


1/2 








(a|+i,+i,+i) = 4 


1 , 
— (ai + 

TT 


102) 


(DIO) 


(a| + i,+i.-i) = - 


TT 




(DU) 


(a| + i,-i,+i) = 4 


i , 
--(flo + 

TT 


103) 


(D12) 


(a|+i,-i,-i) = - 


-(ai - 

TT 


ia2) 


(D13) 



For i 



%/3 



(D14) 



(a|l, +1, +1) = — -j={ax + ia2) 

TTV 2 

/3 

(a|l,+l, 0)=i (ai +m2)(ao + W3) (D15) 

TT 

/3 

(a|l,+l,-l) = ^(ao + ma)^ 



(D16) 



(a|l,-l,+l) 



V3_ 
7r\/2 
.%/3 



(ao - ia3)^ 



(D20) 



(a|l,-l, 0) = i — (fli - m2)(ao - zag) (D21) 



i|l,-l,-l) = 



7r\/2 



(ai - ia2)'' 



(D22) 



For j = 3/2 



V2, 



(D23) 



(a|+i,+i,+i) = [ax+ia^r 

TT 

(a|+|,+i,+|) = i — (fli + ^2)^(00 + ia-i) (D24) 

TT 



2 _i J_3^ 
2 ' 2 2 / 



(fli + ia2)(ao + ia3)^ (D25) 



(a| + |,-|,+|) = -i — (ao + iazY 

TT 



(D26) 



\/6 

= — i — (ai + ia2)^(ao — 103) (D27) 

TT 



\/2 

(fli + ia2)(l - 3(ao + Og)) 



(D28) 



a + 



3 _i ,1 
2 ' 2 2 / 



2 ' 2 2 / 



a/2 

j (ao + ja3)(l — 3(a^ + aa)) 

TT 

(D29) 

— (oi - ia2)(ao + ia3)^ (D30) 

TT 



(a| + |,+f,-i) 



2 ' ' 2 ' 2 , 



— (ai + ia2)(ao - m3)^ 

TT 

.V2, 



(D31) 



= — (ao - ia-i){\ - 3(ai + a^)) 

TT 



(D32) 



3 _i _i 
2 ' 2 ' 2 , 



2 ■ 2 ' 2 / 



= — (ai - m2)(l - 3(a5 + a^)) (D33) 

TT 

= i (ai — m2)^(ao + «a3) (D34) 

TT 
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(D35) 



(a| + f ,+§ -1) = i — {ao - ia,sy^ 

TT 

a/6 

(a|+|,+i,-|) = — (ao - ia3f{ai - ia2) (D36) 



-3 -1 -3 
2 ' 2 ' 2 / 



-I — (ai — 102)^(00 — ias) (D37) 

TT 



/2 

(a|+|,-|,-|) =- — (ai-io2)3 

TT 



(D38) 



those with j = 1, m = 



{R\l, 0,+l) = -(7?3i+«i?32) = -i?°+ (E5) 
{R\l, 0, 0) = i?33 = i?°° (E6) 
(i?|l, 0,-l) = i?3i-ii?32 = i?°- , (E7) 



APPENDIX E: WAVEFUNCTIONS ON THE S'0(3) 
MANIFOLD 



and with j = 1, m 



The lowest state in this case has j = m = n = and 
is simply a constant 



(i?|0,0,0) = 1 . 



(El) 



The next higher states are those with j = 1. There are 
nine of these. Those with j ~ 1, m — 1 

{R\l, +1, +1) = Rii - R22 + i{Ri2 + i?2i) = 

(E2) 

(i?|l,+l, 0) = -(i?i3+ii?23) = --R+° (E3) 
{R\l, +1,-1) = -{Rii + R22 - i{Ri2 - R21)) = -R+- , 



{R\l,-l, +1) = -(i?ii + R22 + i{Ri2 - R21)) = 

(E8) 

(ii|l,-l, 0) = Ri3 - iR23 = R-° (E9) 
{R\l, -1, -1) = Rn - R22 - i{Ri2 + R21) = R— ■ 

(ElO) 



Higher wavefunctions with total (iso)spin j are made up 
of products of i?^^ to power j. The wavefunctions here 



(E4) have not been normalized. 
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